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Abstract 

We analyze the resonant inelastic x-ray scattering (RIXS) spectra at the K edge of Mn in the 
antiferromagnetic insulating manganite LaMnOs. We make use of the Keldysh-type Green-function 
formalism, in which the RIXS intensity is described by a product of an incident-photon-dependent 
factor and a density-density correlation function in the 3d states. We calculate the former factor 
using the 4p density of states given by an ab initio band structure calculation and the latter us- 
ing a multi-orbital tight-binding model. The ground state of the model Hamiltonian is evaluated 
within the Hartree-Fock approximation. Correlation effects are treated within the random phase 
approximation (RPA). We obtain the RIXS intensity in a wide range of energy-loss 2-15 eV. The 
spectral shape is strongly modified by the RPA correlation, showing good agreement with the exper- 
iments. The incident-photon-energy dependence also agrees well with the experiments. The present 
mechanism that the RIXS spectra arise from band-to-band transitions to screen the core-hole po- 
tential is quite different from the orbiton picture previously proposed, enabling a comprehensive 
understanding of the RIXS spectra. 

PACS numbers: 78.70.En 75.47.Lx 71.28.+d 78.20.Bh 
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I. INTRODUCTION 



Recently, the perovskite-type manganites showing the colossal magnetoresistance (CMR) 
effect have attracted much attention because the CMR effect may be the key for high- 
capacity magnetic storage and spintoronics as the next generation electronic devices. It is 
widely recognized that the spin, charge and orbital degrees of freedom of the 3d electrons 
play important roles in the CMR effect in the manganites. In the undoped material LaMnOs, 
the crystal structureii^ belongs to the Pbnm space group below 780K, where each MnOe 
octahedron is noticeably rotated, tilted, and distorted due to the cooperative Jahn-Teller 
distortion (JTD).- The orbital degeneracy in the Cg states is lifted by the JTD or by the or- 
bital exchange interaction similar to the superexchage for spins, ^ forming an orbital-ordered 
state, in addition to an A-type antiferromagnetic long-range order.-"^ With hole doping, a 
variety of spin, charge, and orbital ordered phases appear. It is widely accepted that the 
CMR effect is a consequence of the destruction of the orbital ordered state in the presence of 
a magnetic field, which may break the subtle balance among the multiple degrees of freedom. 

Resonant inelastic x-ray scattering (RIXS) technique using the x-ray tuned at the 
transition-metal K edge has been recognized as a powerful tool to investigate the charge and 
orbital degrees of freedom in transition-metal compounds. The RIXS process is described 
as a second-order optical process, in which a Is-core electron is excited to an empty Ap 
state by absorbing an incident photon, then charge excitations are created in the 3(i states 
to screen the core-hole potential, and finally the photoexcited 4p electron recombines with 
the Is-core hole by emitting a photon. In the final state, the charge excitations are left 
behind in the ?>d states. In contrast to the optical conductivity measurement, it allows us to 
investigate the momentum dependence of the excitations because the corresponding x-rays 
have wavelengths of the same order of lattice spacing. The clear momentum dependence 
have been observed in cuprates.-'^'^'^'ii Note that electron energy loss spectroscopy can also 
detect the momentum dependence, but it cruelly suffers from multiple scattering effects. 

Several RIXS experiments have been carried out on LaMnOs and the hole-doped 
compounds.— "^"^ For the undoped case, the RIXS spectra show three noticeable features 
at the energy loss 2.5, 8, and 11 eV, when the incident photon energy is tuned near the Mn 
K edge.— These features are also observed in the hole-doped manganites.— >^ Based on the 
theoretical analysis exploiting the Liouville operator method,— Inami et al argued that the 
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2.5-eV peak arises from an orbital excitation across the Mott gap, which is driven by the 
offdiagonal part of the Coulomb interaction between a photoexcited 4p electron and an Cg 
electron.— This assumption, however, seems unrealistic, because such a offdiagonal Coulomb 
interaction is much smaller than the Coulomb interaction between the Is hole and 3(i elec- 
trons. This mechanism could not explain the origin of the 8- and 11-eV peaks either. The 
purpose of this paper is to elucidate the origin of the observed features in the RIXS spectra 
of LaMnOs on the basis of the well established assumption that the charge excitations are 
created in the 3(i states to screen the core-hole potential. To obtain the comprehensive 
understanding of the spectra, we use a multi-orbital tight-binding model involving all the 
Mn 3(i orbitals and O 2p orbitals, and take account of the crystal distortion forming Pbnm 
by varying the transfer energy between the 3d and 2p orbitals. 

We calculate the RIXS intensity with the use of the formula proposed by Nomura and 
Igarashi (NI).— The NI formula of the RIXS spectra is a kind of extension of the reso- 
nant Raman theory developed by Nozieres and Abrahams^ on the basis of the many-body 
formalism of Keldysh.— The formula has advantages that it can rather easily be applied to 
complicated models including many orbitals and provides clear physical interpretations of 
the RIXS spectra. This formula utilizes the Born approximation to the core-hole potential, 
and divides the RIXS intensity into two factors; one describes an incident-photon depen- 
dence and the other is the density-density correlation function in the 3d states. Similar 
formulas have been derived by using different methods.— The NI formula has been success- 
fully applied to the quasi-one-dimensional cuprates SrCuOa^^, CuGeOs,^ two-dimensional 
cuprate La2CuO^^'^'^, and the prototypical AFM insulator NiO.— In these studies, the 
calculation have been carried out at zero temperature; the electronic structures in the AFM 
phase have been calculated within the Hartree-Fock approximation (HFA). It is known that 
the HFA works well for the description of electronic structures in the AFM insulators. Two- 
particle correlations have been taken into account within the random phase approximation 
(RPA). It has been found that the RPA correction modifies strongly the spectral shape as 
a function of energy loss, having led to a good agreement with the experiments. With these 
successes, we may conclude that the RIXS intensity arises from band-to-band transitions to 
screen the core-hole potential in the intermediate state. Multiple-scattering contributions 
due to the core-hole potential have been also investigated in order to examine the validity 
of the Born approximation, because the core-hole potential is not definitely weak.— Having 
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evaluated the contributions by means of the time-representation method by Nozieres and 
De Dominicis,— it was found that the contributions could be mainly absorbed into the shift 
of the core-level energy with minor modifications of the RIXS spectral shape.— This result 
partly justifies the use of the Born approximation. 

In the present study of LaMnOs, we treat the strong Coulomb interaction between the 
3d orbitals within the HFA; we obtain an A-type AFM insulating solution with an energy 
gap ~ 1.0 eV, where the occupation on the 3z^ — r^-type orbital is larger than that on the 

— y^-type orbital at each Mn site. This result corresponds well to the observed orbital 
ordered state. Note that the band structure calculation with the local density approximation 
(LDA) fails to reproduce the wide energy gap. We calculate the density-density correlation 
function at zero temperature by using the energy bands thus obtained and by treating the 
two-particle correlations within the RPA. We calculate another factor, the incident-photon- 
dependent factor, using the 4p density of states (DOS) obtained from the ab-initio band 
structure calculation. Combining the two factors, we finally obtain the RIXS spectra, which 
show good agreement with the experiments.—'^'^ We could assign the 2.5-eV peak in the 
RIXS spectra as a function of energy loss to the electron-hole excitation across the Mott gap 
in local majority spin channel, and the 8-(ll-)eV peak to the transition from the occupied Cg 
states strongly hybridized with the O 2p states in the deep valence band to the unoccupied 
Cg states in the local majority (minority) spin channel. We also make clear the origin of 
the incident-photon-dependence of the spectra by examining the corresponding factor in our 
formula. 

The present paper is organized as follows. In Sec. Ill Al we introduce the multiorbital 
tight-binding model. In Sec. Ill Bl we discuss the electronic structure within the HFA in 
the AFM phase of LaMnOs. In Sec. IIIII we briefiy summarize the NI formula for the 
RIXS spectra. In Sec. HVl we present the calculated RIXS spectra in comparison with the 
experiments. The last section is devoted to the concluding remarks. 
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Figure 1: (a) Sketch of the unit cell. La sites are omitted. The experimentally observed crystal and 
magnetic structure is assumed. The arrows on the Mn sites indicate the magnetic moment, (b) 
Local coordinates at each Mn site. The z\ (A = 1, • • • ,4) axis is taken to be parallel to the longest 
Mn;s^-0 bond, while x\ and y\ axes are nearly parallel to the shortest Mn;^-© and the middle length 
Mn;s^-0' bonds, respectively. 



Table L Tight-binding parameters in units of eV. SK parameters for the nearest Mn-0' pair and for 
the nearest 0-0' pair are given in the table, while the other parameters are determined by adopting 



° law suggested by Harrison.— Slater Integral and are taken from ref. 
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5.03 


U 4.5 


pdn 


0.93 




9.74 




ppa 


0.63 




6.96 




ppn 


-0.18 


Ed — Ep 


-14.0 


A 4.0 



We assume the crystal and magnetic structure of LaMnOa as shown in Fig. [T]-a and also 
define the local coordinates x\, y\, and zx on the A-th Mn site in a unit cell as shown in fig. 
[l]-b. We introduce a tight-binding model involving all Mn 3(i orbitals and O 2p orbitals. 



We exclude orbitals belonging to the La atoms, since those orbitals play minor roles in the 
electronic states near the insulating gap. Thus, the model Hamiltonian is expressed as 

H = Ho + Hi, (1) 



ima jma 

+ E E {^tm,jm'4maPjm'a + H.C.) 
(i,j) mm'cr 

+ E E {^7rn,l'rn'P]maPj'm'a + H.C.) , (2) 
mm' a 

= ^ E E ^ '^^i' ^^3, ^i) dl^dl/^J^,. (3) 

The part Hq is the kinetic energy part, where dima and Pjma denote the annihilation op- 
erators of an electron with spin a in the 3d orbital m at Mn site i and in the 2p orbital 
m at O site j, respectively, nf^^ and n^^^ are the number operators given by d\^^dima 
and p]rnaPjrnc7-, respectively. The transfer integrals, f^^ -,^, and t^^j,^, are evaluated from 
the Slater- Koster (SK) two-center integrals, (pda), (pdn), (ppcr), (ppyr).— We neglect the 
hybridization between the Mn sites. In order to take account of the JTD, we assume the 
law for the SK parameters suggested by Harrison,— where i represents the atomic distance 
between neighboring sites and a = 2 for ppa and ppvr, a = 3.5 for pda and pdn. Thus, the 
JTD, rotation, and tilt of the MnOg octahedron are incorporated into the model. The part 
Hi represents the intra-atomic Coulomb interaction on the Mn sites. The interaction matrix 
element g (z/i, z/2; 1^3, 1^4,), where u stands for spin-orbit (ma), is written in terms of the Slater 
integrals F^, and F*. Among them, and F'^, which are known to be slightly screened 
by solid-state effects, are taken from the cluster model analysis of the x-ray photoemission 
spectroscopy.-2S On the other hand, F° is known to be considerably screened, so that we re- 
gard the value as an adjustable parameter. The Coulomb interaction on O sites is absorbed 
into a renormalization of the O 2p level parameters Ej^^. The Mn d-level position relative 
to the O p-levels is given by the charge-transfer energy A defined as A = Ed — Ep + AU 
in the d^ configuration, where U is the multiplet-averaged d-d Coulomb interaction given 
by f/ = — (2/63) — (2/63) F"^.— The charge transfer energy A is also treated as an 
adjustable parameter in our calculation. The parameters used in the calculation are listed 
in Table [D 
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B. Hartree-Fock approximation 

Assuming the A-type AFM order, we solve the tight-binding Schrodinger equation within 
the HFA. We obtain a stable A-type AFM solution, which has the energy gap ~ 1.0 eV and 
the spin moment at Mn site ~ 2.0 h. The densities of states (DOS) projected onto the 
3zl — r^, x\ — yl, and t2g states are shown in fig. [2], and that projected onto 2p states in fig. 
[HI In the local coordinates, they are independent of A (1, ■ ■ ■ ,4). The local majority spin 
3zl — states are almost fully occupied and mainly concentrate on the occupied energy 
ranges denoted by A and C; they also concentrate on the unoccupied energy range denoted 
by A'. The local majority spin x\ — y\ states are partially occupied on energy ranges 
denoted by A and C; they also concentrate on the unoccupied energy ranges denoted by A'. 
On the other hand, the local minority spin states highly concentrate on the unoccupied 
energy range denoted by B' and C, although small amount of them resides on the occupied 
energy range denoted by B. The difference in the occupation numbers between the ?>z\ — 
and x\ — y\ states implies that the eg states are orbitally ordered, corresponding to the 
experimentally observed orbital order. In contrast to the Cg states, the local majority and 
minority spin t2g states are almost perfectly occupied and unoccupied, respectively. The Cg 
states around the energy ranges denoted by B and C are highly hybridized with O 2p states. 
In contrast, the occupied states around the energy range denoted by A and unoccupied states 
denoted by A' have very small weight of the 2p states. This difference arises from the fact 
that the oxygen 2p states are located in the relatively deep energy region as shown in fig. [3l 
LaMnOs is close not to the charge-transfer-type insulator but to the the Mott-Hubbard-type 
insulator. 

The dispersion curves along some symmetric lines near the gap are shown in fig. [H Labels 
assigned to the curves correspond to the states shown in the figO The curves corresponding 
the states denoted by A and A' well reproduce those calculated with the ab initio band 
structure calculation based on the LDA+[/ method^ except for the magnitude of the gap. 
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Energy (eV) 

Figure 2: DOS projected on the Mn — r^, x\ — and t2g states. The origin of energy is at the 
top of valence band. In each panel the upper half represents the DOS for the local majority spin 
and the lower one for the local minority spin. 
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Figure 3: DOS projected on the O 2p states at the sites denoted by O and O' in fig. [H The origin 
of energy is at the top of valence band. 
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Figure 4: Dispersion curves along the several symmetric lines vicinity of the gap. The origin of 
energy is at the top of valence band. 



III. RIXS PROCESS 



We briefly summarize the NI formula for the RIXS, following ref. l23|. In the RIXS process, 
the incident photon is absorbed by exciting a Mn Is core electron to the unoccupied Mn 
4p state, and a photon is emitted by recombining the 4p electron and the core hole. This 
process may be described by 



w 



(4) 



where w represents the dipole transition matrix element between the Is and the 4p states. 
We assume that w is constant, since it is expected to change little in the energy range 
of 20 eV above the absorption edge. It is known that the energy dependence of w gives 
minor effect on the Cu i^-edge absorption spectra.— The represents the 77-th component 
{rj = X, y, z) of the photon polarization vector. Annihilation operators p'^^^ and Sia- are for 
states Apr, and state Is with spin a at Mn site i , respectively. The annihilation operator 
is for photon with momentum q and polarization In the intermediate state of the 



RIXS process, the core-hole potential is acting on the 3(i states, creating an electron-hole 
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Figure 5: Diagram for the RIXS intensity within the Born approximation for the Is core-hole poten- 
tial. The wavy and dotted lines represent photon Green's functions and the core-hole interaction 
V, respectively. The solid lines with the labels Ap and Is represent the bare Green's functions 
for the Ap electron and the Is core electron, respectively. The elliptic part between the core-hole 
interaction lines corresponds to the density-density correlation function of the Keldysh type. The 
lines with double arrows are the Keldysh type Green's functions. The shaded area represents the 
effective scattering vertex renormalized by 3d-3d Coulomb interaction in the RPA. 

pair within the Born approximation. The interaction is described as 

Hls-3d = y 4madimas\^,Siai, (5) 

imaa' 

where i runs over Mn sites. Note that, although the value of the core-hole interaction V is 
not known and may strongly depend on the model, it is expected to be much larger than the 
4p-3(i offdiagonal Coulomb interaction. In the end of the process, an electron-hole pair is 
left behind carrying momentum-energy q = (q, cu) = (q^ — q/,^^^ — ujf), where qi = {qi,uJi) 
and qj = (q/, ujf) are momentum-energy of incident and scattered photons, respectively. 

The RIXS intensity is derived on the basis of the Keldysh-Green function scheme. It is 
diagrammatically represented in FiglHl Within the Born approximation for the core-hole 
potential, we obtain 

W{q., e,,; qj, ej) = Yl E ^^''■^"'""^'^nt7.,A',nV'(?)'/BAA' (^., e.; ujf, e;) , (6) 

\ma X'm'a' 

where A indicates the A-th Mn site in a unit cell, ua represents the position vectors of the 
A-th Mn site in a unit cell. N is the number of unit cells. The factor Jbaa' ("^i, Gj; cj/, e/) 
describes the incident-photon dependence, which is given by 

Jbaa' (t^i, Gf, ujf,ef) = I ^ e^^L"^^ i^^u ^/) e/v j ^^^v^bx' i^i'^ ^f) ^fv'\ ' C^) 

\ qri' / \ rjrj' J 
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where e,^ (e/r?) is the 77-th component of the polarization vector (ej) with rj = x,y, z, and 
Lbx i^h^) is given by 



The Tis represents the life-time broadening width of the core-hole state, and the lower limit 
of the integral eo indicates the energy at the bottom of the 4p band. This expression comes 
from the upper triangle in FiglSl The pj^^ is the DOS matrix in the p symmetric states at 
the A-th Mn site, which may be given by 

Pit (e) = E E <^A,. in, k) in, k) <5 (e - (k)) , (9) 

(T n.k 

where (pxr^a in, k) is the amplitude of Pr, component with spin a at the A-th Mn site in the 
band state specified by the band index n and crystal momentum k with energy e„ (k). The 
factor Yy^,^, (g) in Eq. (I6|) is the density-density correlation function of the Keldysh 
type, which is defined by 

/oo 
((PqA',n'.')^(^)PqA™.(0))e*"Mr, (10) 
-oo 



where 



with 



PciXma — \/ T7 E'^k+qAma'^kAmCT, (H) 



d 



kAmcr 



The index \ma specifies a tight-binding orbital at site A with orbital m and spin a; A 
(= 1,2,3,4) is assigned to four Mn sites in a unit cell. Wavevector k in eq. (fTTl) runs over 
the first Brillouin zone. Vector r„ in eq. (fT2ll represents a position vector of the n-th unit 
cell. A single phase factor k ■ r„ is assigned to all the 3d states in each unit cell, and hence 
an extra factor e*'!'^"^""^'^ is required in eq. (jGj). 

We calculate the density-density correlation function ( ITOl l at zero temperature with taking 
account of the correlation effect on the electron-hole pair by the RPA. Abbreviating the 
indices \ma as ^, it may be written as 

(^) = E E ^hu' (^) KiL (^) ii) , (13) 
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where 



n 



-(0) 



UJ 



E,, (k + q) + E, (k)) [1 - n.y (k + q)] n, (k) 



V'Sij' (k + q) (k + q) </?5^j (k) (k) . 



(14) 



Ej (k) and (k) are the eigen-energy and the occupation number of the eigenstate specified 
by jk, respectively. The (pxma,j (k) represents the amplitude of the 3d orbital and spin ma 
at the Mn site A in the energy eigenstate specified by jk within the HFA. The RPA vertex 
is given by 



(15) 



where / represents a unit matrix, and f is the bare four-point antisymmetric vertex given 

by 



The two-particle propagator F (g) is given by 

1 

-I 5l52iS354 



{g (z/iZ/2; Z/3Z/4) - g (z/iz/g; i^ii^3)) Sx^xJxaxJxiXs- 



(16) 



^ E E ^^€4 J (k) (k) <^5,,,' (k + q) ifilj, (k + q) 



X 



rij (k) [1 - nj> (k + q)] 



Uji (k + q) [1 - (k)] 



(17) 



u - Ej, (k + q) + Ej (k) +i5 uj - Ej, (k + q) + Ej (k) - i5 _ 
We evaluate the density-density correlation function f fTOl) using eqs. (fT3l)-(fT7l). For more 



details of the derivation, see refs 



21 



23l . and 



22 



Note that the terms involving n~+*^°) (g) are 
neglected; they have no contribution for > at zero temperature, because n~+*^°) (g) oc 
X^k (-'- — (k)) n (k + q) (5 (cj — Ej, (k + q) + Ej (k)) with omitting unimportant factors. 



IV. CALCULATED RESULTS 

In order to calculate the incident-photon-dependent factor Jbaa' i^i^^u^f^^f), we need 
the Ap DOS Pxipi^) on the A-th Mn site. We evaluate the Ap DOS using full-potential-linear- 
augmented-plane-wave band structure calculation based on the LDA+f/ method. It should 
be noted here that the Mn 4p DOS depends on r], rj' and on Mn site A due to the strong JTD, 
giving rise to the Mn i^-edge resonant elastic x-ray scattering intensity on forbidden Bragg 
spots.— "^"^ Figure [6] shows the 4p DOS averaged with t] {J2ri Pxipi^)) convoluted with a 
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Lorentzian function of FWHM 2Tis = 2eV, in comparison with the absorption experiment.— 
We set the energy difference between the Mn Is level and the prominent peak in the 4p DOS 
to be 6554 eV. Under the condition that the dipole matrix element is constant and that the 
interaction is neglected between the core hole and the 4p electron, the 4p DOS becomes 
proportional to the Mn i^'-edge absorption spectra. The agreement with the experiment 
indicates that the above condition is nearly satisfied. 

Equations ((6])- (El) indicate that the photon polarization dependence of the RIXS intensity 
correlates only with the Mn 4p states. Therefore, the polarization dependence would not 
have relevant information on the charge excitations in the 3d states. Hence, we will not touch 
on the polarization dependence. Assuming polarization-unresolved analysis, we neglect the 
offdiagonal elements of the 4p DOS matrix and replace each diagonal element with the 
averaged value for simplicity. Accordingly, the factor Jbaa' (t^i, e^; cu/, ej) is replaced with 

2 

the averaged site-independent factor JxioJi^Uf) ~ X]?, -^ba (^«' ^/) ■ Figure [7] shows the 
contour plot of Ja {uJi]UJf) as a function of incident photon energy uji and energy loss u = 
uji — ujf. The enhancement appears around the energy loss of 2 eV for the incident photon 
energy around 6554 eV. As the incident photon energy increases, the enhancement peak 
moves toward higher energy loss region with decreasing the intensity. It is expected that 
the factor Jbaa' (t^i, cu/, ej) has dependence similar to the factor Ja (cu,; cu/), although 
the enhancement peak position and the intensity may somewhat depend on the photon 
polarization. 

Another factor Yxim' a XmS'l) mainly determines the structure of the RIXS spectra as a 
function of energy loss. We calculate this factor from eqs. (|T0l)-(fT7ll using the tight-binding 
wavefunction within the HFA and by including the RPA correction. The transition from the 
occupied Cg states in the energy ranges denoted by A, B and C to the unoccupied Cg states 
in those denoted by A', B', and C can be driven due to the Is core-hole potential, because 
the occupied states and the unoccupied states comprise the Cg states with the same local 
symmetry. Therefore, the Cg states can contribute to the RIXS intensity. In contrast to the 
Cg states, the t2g states cannot contribute to the RIXS intensity, since the local majority and 
minority spin t2g states are almost perfectly occupied and unoccupied, respectively, and the 
spin-flip transition is not allowed. 

Combining two factors, we obtain the RIXS spectra. We convolute the result with a 
Lorentzian function of FWHM = 0.5 eV for taking account of the instrumental resolution. 
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Figure 6: The LDA 4p DOS convoluted by the Lorentzian function with FWHM = 2eV (thick 
sohd hne curve). In connection with the absorption coefficient, the origin of energy is shifted so 
that the prominent peak locates at 6554 eV. Thin line curve is the absorption spectra reproduced 



from ref. 
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Figure 7: Contour plot of Ja (ajj; w/) as a function of energy loss uj = uJi — ujf and incident photon 
energy Wj. 2ris = 2 eV. Units of the intensity are arbitrary. 

Figure [8] shows the spectra thus evaluated at the momentum transfer q = (1.6,1.6,0.0) 
as a function of energy loss for several incident photon energies in comparison vv^ith the 



experiments.- 



IJ.I:;.I4 



We obtain continuous spectra ranging from u; = 2 eV to 15 eV vv^ithin the 



HFA; there are three features around u = 2.5, 7.5, and 12 eV, vi^hich arise from transitions 
of A — i> A', C ^ A' in the local majority spin channel, and B B', C in the local minority 
spin channel, respectively, vv^ithin the Cg states (see Figs. [2]and[ll). As already stated, the 
t2g states could not contribute to the RIXS intensity. The spectral shape vv^ithin the HFA 
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is drastically modified by the RPA correction; the intensity around the energy loss uj ~ 2.5 
and 8 eV is suppressed and that around ~ 12 eV is enhanced. The spectral shape thus 
modified corresponds well to the experimentally observed features at = 2.5, 8, and 11 eV. 
Another characteristic is that the weight of the RIXS intensity moves toward higher energy 
loss region as the incident photon energy increases. This change is nicely reproduced by the 
calculation, mainly due to the effect of the incident-photon-energy factor Ja. Figure [9] shows 
another comparison with the experiment at other momentum transfers q = (2.7, 0.0, 0.0) and 
(4.5,0.0,0.0).— The calculated spectra agree with the experiment. 
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Figure 8: RIXS spectra at q = (1.6, 1.6, 0.0) as a function of energy loss for several incident photon 
energies. Calculated RIXS spectra are shown in the left panel and the observed spectra reproduced 



from ref. |U in the right panel. The solid and dashed curves are spectra calculated by including the 
RPA correction and within the HFA, respectively. The incident photon energies are shown besides 
each spectral curve in units of eV. Arrows indicate characteristic features in the spectra. 
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Figure 9: RIXS spectra at q = (2.7,0.0,0.0) and (4.5,0.0,0.0) as a function of energy loss for 
several incident photon energies. The solid curves represent the spectra calculated by includin g th e 



RPA correction. Dotted ones are the observed spectra at hwi = 6556 eV reproduced from ref. 
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FigurefTOlshows a contour plot of RIXS spectra as a function of energy loss and momentum 
transfer q along (0, 0, 0)-(2, 0, 0), -(0, 2, 0), -(0, 0, 2), and -(2, 2, 0) lines. The incident photon 
energy cUj is fixed at 6556 eV. Three ridges around the energy loss 2.5, 7.5, and 12 eV 
correspond to the three peaks around the energy loss 2.5, 7.5, and 12 eV, respectively, in 
figs. [8]and[9l The peak position only slightly depends on the momentum transfer q, moving 
within ~ 1 eV. This is consistent with the experiments and may be due to the fact that 
the states relevant to the RIXS process concentrate on the narrow energy ranges. Note 
that the peaks at 2.5, 7.5, and 12 eV correspond to transitions A ^ A', C — A' in the 
local majority spin channel, and B B', C in the local minority spin channel, respectively 
(see Figs. [2]and[l]). In contrast, the peak intensity noticeably depends on the momentum 
transfer showing the period of 2 along the a, b, and c directions, although this fact has 
not been emphasized in the experiments.—'^'^ It does not seem easy, however, to explain 
the origin of this dependence, since the intensities depend sensitively on the 3d weights in 
the energy bands (see Eq. (fT4l)). Different from the present case, the peak positions as a 
function of energy loss have been found clearly moving with varying momentum transfer in 
La2Cu04 in the high-resolution experiments.-'^'^'^'^ It should be noted here that the peak 
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shift with varying momentum transfer should not be interpreted as a dispersion relation 
of a kind of exciton but as a change of spectral weight in the continuum spectra. The 
crystal momentum dependence of the weight of the states relevant to the RIXS process 
and the energy dispersion of the single electron states determine the momentum transfer 
dependence of the RIXS intensity. 




Figure 10: RIXS spectra as a function of momentum transfer q along (0, 0, 0)-(2, 0, 0), -(0,2,0), 
-(0,0,2), and -(2,2,0) lines. Pbnm setting is assumed for the definition of q. The incident-photon 
energy is fixed at = 6556 eV. The unit of the intensity is arbitrary. 



V. CONCLUDING REMARKS 

We have analyzed the incident-photon-energy and momentum dependence on the RIXS 
spectra as a function of energy loss in LaMnOs. We have utilized the formula developed by 
Nomura and Igarashi, which expresses the RIXS spectra by a product of the photon depen- 
dent factor and the density-density correlation function. The former factor, which describes 
the dependence on the incident-photon energy and polarization, has been calculated by us- 
ing the Mn 4p DOS given by the ab initio band structure calculation. The latter function 
describes the charge excitations generated in the 3d states to screen the core-hole potential 
in the intermediate state. 

Having combined the above factors, we have calculated the RIXS spectra as a function 
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of energy loss. We have obtained a continuous spectrum consisting of three features in good 
agreement with the experiments. We have found that the RPA correlation modifies drasti- 
cally the spectral shape, which indicates the importance of electron correlations for making 
quantitative analysis. We have demonstrated that the peak intensity depends strongly on 
the momentum transfer q with the periodicity 2 along the a, b, and c directions while the 
three peak positions as a function of energy loss only slightly depend on q. Furthermore, 
we have obtained the spectral weight moving toward high energy-loss region with increasing 
incident-photon energy, in agreement with the experiments. 

The present analysis naturally leads to a RIXS picture that the spectra are brought about 
by band-to-band transitions (augmented by the RPA correction) in order to screen the core- 
hole potential. This picture is different from a previous picture that the the spectra are 
brought about by 3d orbital excitations (orbitons) created with the use of the offdiagonal 
part of the Coulomb interaction between the photo-excited 4p electron and 3d electrons.—'^ 
This "orbiton" picture seems unreasonable, since the 4:p-3d offdiagonal Coulomb interaction, 
which causes orbitons, is much smaller than the ls-3d Coulomb interaction, which causes 
band-to-band transitions. Actually the "orbiton" picture was not successful in providing the 
spectra comparable to the experiment.— Note that the two pictures give different selection 
rules for the creation of excitation; the present picture forbids the transition from 3z^ — re- 
type orbital to the — y^-type orbital, while the latter allows it by changing 4p states. 

For doped cuprates and manganites, experimental data have been accumulated.—'^ In 
doped cuprates the RIXS spectra have been analyzec _ 
within the HFA.— As already demonstrated in Refs. 
taking account of the RPA correlation, is not sufficient for quantitative understanding of 
the spectra even in the undoped materials. In doped materials, electron correlations may 
become more important, and the HFA-RPA scheme would not work well. Also, the Born 
approximation may be insufficient for treating the core-hole potential, since many electron- 
hole pairs could be created in the absence of the energy gap. It seems hard to answer these 
questions by analyses with a detailed model like the present paper, and such studies are left 
in future. 
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